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Abstract. We show that, for the Id cubic NLS equation, widely separated equal 
amplitude in-phase solitons attract and opposite-phase solitons repel. Our result 
gives an exact description of the evolution of the two solitons valid until the soli- 
tons have moved a distance comparable to the logarithm of the initial separation. 
Our method does not use the inverse scattering theory and should be applicable 
to nonintegrable equations with local nonlinearities that support solitons with ex- 
ponentially decaying tails. The result is presented as a special case of a general 
framework which also addresses, for example, the dynamics of single solitons sub- 
ject to external forces as in [3]. 



1. Introduction 
We consider the Id nonlinear Schrodinger equation (NLS) 
(1.1) id t u + \d 2 x u+ \u\ 2 u = 0. 

It has a single soliton solution u(x,t) = e lt ' 2 sechx. The invariances of ( |1.1[ ) can be 
applied to produce a whole family of solutions. To describe them, let 

^^"^-"Vsech^z-o)) 



(1.2) 



rj(x, fi, a, 0, v) 



e e ■ 



for parameters 0, a, v G K, fi > 0. Then u(x, t) = r/(x, fi(t), a(t), 0(t), v(t)) solves (1.1 ) 
provided 

a(t) = a + tv fj,Q l 



(1.3) 



K V(t) = Vq 



In this paper, we study the evolution of initial data that is the sum of two widely 
separated solitons: 



(1.4) 



u Q (x) = T](x, fl 10 , a w , 0io, Vxo) + T)(x, /i 2 0, O 2 0, 6*20, v 2 o ) 



1 We order the parameters as (ji,a,9,v) to mimic (q 1 , Q 2 , p 1 , p 2 ) as canonical coordinates for the 
four dimensional symplectic space with symplectic form dp 1 A dq 1 + dp 2 A dq 2 . 



(1.5) tto(ar) 
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where \a 20 — «io| ^ !• I n particular, we focus on two illustrative cases. In both cases, 
we consider identical mass solitons with zero initial velocity. In Case 0, we take the 
same initial phase, corresponding to an even superposition and in Case 1, we take 
opposite initial phase corresponding to an odd superposition. 

rj(x, 1, — do, 0, 0) + r)(x, 1, ao, 0, 0) Case a = 
r](x, 1, — a , re, 0) + r](x, 1, a , 0, 0) Case a = 1 

We find that in the same phase case (Case 0), the two solitons are drawn toward 
each other and in the opposite phase case (Case 1) they are pushed apart- see Fig. 
In either case, the solution u to (1.1) is well- approximated by 

(1.6) u z {x) = 7](x, n x , ax,9i,vx) + r](x, fi 2 , a 2 , 9 2 , v 2 ) 
where z represents coordinate^] 

(1.7) z = (z 1 , z 2 , z 3 , z 4 , z 5 , z 6 , z 7 , z s ) = (fix, ax, /x 2 , a 2 , Ox, vx,9 2 , v 2 ) 



As parity is preserved by the flow (1.1), we have 

/., „\ def def def 

(1.8) fl = fix = A*2 , a = — Ol =0,2, V = —V\ = V 2 , 

and 9 = f 9x = 2 in the same phase case (Case 0), while 9 = f Ox — n = 2 in the 
opposite phase case (Case 1). 



Theorem 1.1. Suppose that u(t) is the solution to (1.1) with initial data (1.5). Let 
h = e~ a ° < 1 (so a = \ogh~ 1 > I). Let 




Case a = 
Case a = 1 




T 

Let (a{t),v{t)) solve 
(1-9) 

1 • = -4{-iy e 
with initial data (a ,0). Let ji solve 

(1.10) fi= (-If (&a - A)ve~ 2a , 
and then let 6 solve 

(1.11) = + + 18(-l) ff e- 2a 
Then on <t <T, we have 

\\u(t) -u z \\m < h 2 - , 



2 Supcrscripts are used on z to conform with geometric summation conventions used later in the 



paper. 
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Figure 1.1. The top plot is a depiction of Case (same phase; even 
solution), where the two solitons are pulled toward each other. The 
bottom plot depicts Case 1 (opposite phase; odd solution), where they 
repel. In each case, the solution is modeled in Theorem as u ~ u z = 



r](fi, —a, 9 + an, —v) +r)(fj,, a, 9, v) where (/i, a, 9, v) solve a specific ODE 
system. 



where 
(1.12) 



u z = r](fi, —a, 9 + an, —v) + r](fi, a, 9, v) 



Let us make some remarks on the ODE system (1.9). The energy associated to 
this system is 

H cS = v 2 - A(-l) u e- 2a = -A(-\ye- 2a ° 
In the case a = (same phase), we have v < 0, a < a and 

a(t) = cto — logsec(2/it) 
v(t) = 2/itan(2M) 

valid for < t < e a ° = hr x . In the case a = 1 (opposite phase), we have v > 0, 
a > a and 

a(t) = a + log cosh(2/it) 
v(t) = 2htaah(2ht) 
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valid for < t < a e a ° = ft. -1 log /i -1 . In either case, \x evolving according to (1.10) 
satisfies \\l — 1| < /i 2 and can thus be replaced by 1 in ( |1.12[ ). However, the 0(h 2 ) 
behavior of fi is dynamically significant in that it yields 0(h l ) effects in 9 through 



(1.11). It is evident from the explicit forms for a(t) given above that, on the indicated 



time scale, the soliton has moved a distance comparable to loga . 



We remark that although (1.1) is completely integrable, we do not use the inverse 



scattering theory of Zakharov-Shabat |20j. We expect that one could compute the 
scattering data associated to our initial condition and conduct an analysis using in- 
verse scattering theory that would describe the dynamics for all time. Our argument, 
however, has the merit of being relatively simple and should adapt to most nonin- 
tegrable nonlinearities that support stable solitons with exponentially decaying tails. 
An important example of such a nonintegrable equation is the Id cubic-quintic NLS: 



1 

2 ( 



id t u + -d\u + \u\ 2 u — e\u\ u 







Furthermore, our goal was not just to obtain Theorem |1.1| but to present it in the 
conceptual (yet rigorous) framework of symplectic restriction that illustrates its con- 
nection to previous work of the first author, Holmer-Zworski [7JIE]- 

We cite two papers from the physics literature as motivation for our problem. 
Stegeman-Segev [17] provide an overview of phase-driven two-soliton interaction in 



the context of optics, beginning with an account of the Id case (1.1) that we study 
(see their Fig. 4) and proceeding to a discussion of two-soliton interaction in two 
dimensions in which the attractive forces between in-phase solitons can lead to spi- 
raling structures - see their Fig. 6. The NLS equation also arises in a completely 
different physical setting, Bose-Einstein condensation. Strecker et.al. [IB] describe an 



experiment producing multiple solitons, in which the model is (1.1) with a confining 



potential. A train of five solitons with successively opposite phases are produced and 
oscillate in a well. At the peak of the oscillations, the solitons bunch up but retain 
some separation; [18] explains this in terms of their phase differences. 



We will now give an explanation of Theorem 1.1 and an overview of the proof. 
Consider L 2 (R; C) manifold with metric 



9u(vi,v 2 ) 



(vi,v 2 ) = Re 



v 1 v 2 for u G L 2 , Vi,v 2 G T U L 2 



Introduce J — —i, viewed as an operator T U L 2 — y T U L 2 . The corresponding symplec- 
tic form is 



(1.13) 



u u (vi,v 2 ) = (vi, J 1 v 2 ) 
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Take as Hamiltonian the (densely denned, with domain^] D = H 2 ) function H : L 2 
K given by 

(1.14) = \J \u x \ 2 -\J H 4 

Then 

#'(u) G T U *L 2 ~ T U L 2 

metric g 



The corresponding flow is dtu = JH'iu) yielding (1.1). 



Recalling that rj is given by (1.2), consider the manifold of solitons 
M = { rj(; fi, a, 6, v) | fi > 0, 9 E R, a E R, v G R } . 
Computations show that the restriction of the symplectic form u to M is 

i*u = d6 A dfi + dv A da , 
while the restriction of the Hamiltonian H to M is 



6 



^ 3 , 



Note that the free single soliton flow (1.3) is just the solution to the Hamilton equa- 



tions of motion for H{rj) with respect to i*uj: 

( fi = d e H( V ) = 
a = d v H(rj) = fi~~ 1 v 

6 = -d,H( V ) = l -^ 2 v 2 + 
, v = -d a H{r]) = 

Turning to the double soliton problem, recall that we model the u in terms of u z 



given by (1.6) where z 
notation 



(z 1 , . . . , z 8 ) is given by (1.7). We introduce the shorthand 

def 



Vj = v(-i to, a h Oji v j) > 3 = h 2 - 



Also recall that h = e "° ■€ 1, and the initial soliton separation is 2ao = 2 log h 1 ^> 1. 
Expanding the nonlinearity, we obtain 

(1.15) H p (u z ) = H p ( Vl ) + (H' p (rj 1 ),rj 2 )+H p ( V2 ) + (H' p (rj 2 ) , Vl ) + 0(h 4 ) 

The last two terms are dominant near a 2 (on the effective support of r} 2 ), so that the 
second soliton sees an "effective" Hamiltonian 

(1.16) H eS (fx 2 ,a 2 ,9 2 ,v 2 ) = H( m ) + (H' p (rj 2 ) , Vl ) 

3 This domain is chosen so that JH'(u) — —i{—\u xx — \u\ 2 u) g L 2 . Although we restrict to 
u € D = H 2 here, we will prove estimates on the corresponding flow in H 1 . This parallels the 
situation in the theory of linear self-adjoint operators A, where a dense domain is specified but the 
flow associated to —iA extends to a unitary operator on all of L 2 . 
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and thus its expected equations of motion are 

' fi 2 = do 2 H(r) 2 ) + d d2 (H' p (r] 2 ),r] 1 ) 
a 2 = d V2 H(r) 2 ) + d V2 (H p {r) 2 ),r)x) 
6 2 = -d, 2 H(r ]2 )-d, 2 (H' p (r ]2 ),r ]l ) 
{V2 = -d a2 H(r] 2 ) - d a2 (H' p {r] 2 ), rjx) 



(1.17) 



Likewise, the first two terms in (1.15) are dominant near a\ so the first soliton sees 
an effective Hamiltonian 

H eS (fi 1 ,a 1 ,6 1 ,v 1 ) = H (Vi) + (H p (r]i),m) 

and thus its expected equations of motion are 

* fix = de 1 H(r) 1 ) + d 6l (H p (r]x), 772) 

di = d Vl H(rjx) + d vl (H p (rix),r} 2 ) 

9 1 = -d, 1 H( Vl )-d IM1 (H p (rj 1 ),r l2 ) 

L *i = -daiHirjx) - d ai (H'(rii),ri2) 



fl.181 



Pulling (1.17) and (1.18) together gives us a systems of eight equations in eight 
unknowns: 

' fa = de^H^x) + d dl {H' p {Vi)^2) 
di = d^H^x) + d vl (H p (r]x),r} 2 ) 

e 1 = -d, 1 H( Vl )-d, 1 (K(vi),V2) 

vx = -d ai H{r)x) - d ai (H' p (r]x),m) 
fi 2 = de 2 H(r] 2 ) + do 2 (H p (r} 2 ), rji) 
a 2 = d V2 H(r) 2 ) + d V2 (H p (r] 2 ) , rjx) 
9 2 = -d, 2 H(r ]2 )-d, 2 (H' p ( V2 ),r ]l ) 
v 2 = -d a2 H(r] 2 ) - d a2 (H' p (r) 2 ),r]x) 



(1.19) 



After the even/odd symmetry assumption is imposed, one can distill from (1.19) the 
equations appearing in the statement of Theorem |l.l[ ^ 

We find that the above argument yielding ( 1.19[ ) is a little too vague to adapt to a 
rigorous proof. We now consider a different perspective that informally produces the 



In fact, the above heuristic argument does not invoke the even/odd symmetry assumption and 
thus we might expect the equations (1.19) even without this assumption. However, the equations 
( 1.19 ) are only expected to be accurate to order 0(h 4 ). In the presence of the symmetry assumption 
the eight equations in (1.19) dramatically decouple as (1.9), ( 1.10[ ), ( | 1 . 1 1 [ ) which permits a direct 
analysis of these ODEs that shows that an 0(h 4 ) unknown can only only have a limited 0(h 2 ) effect 
the solution. In the general case, the eight equations in (1.19) are more interdependent and we are 
not certain as to the effect of 0(h 4 ) perturbations. This is not the only obstacle to removing the 
symmetry assumption; see comments below. 
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same set of equations (1.19) but adapts to yield a proof of Theorem and in fact 
extends and unifies the results of El El E] • Recalling z defined in ( |1.7[ ), consider 
now the eight-dimensional two-soliton manifold 

M = { U z = T]i + 7] 2 = T}(; /ii, X , 6>i, V X ) + 7/(-, /i 2 , 2 , #2, ^2) } 



The symplectic form (1.13) restricted to M is 

1 8 

(1.20) i* u = - 2^ ^mO) ^ A cfc r 

where 



'm=l 



= (a^ m (2;)) , a em (z) = (d z m z , J l d z mU z ) 



Let H{u z ) denote the restriction to M of the Hamiltonian (1.14). The expected 



equations of motion for z m are Hamilton's equations for H{u z ) with respect to i*u. 
These equations are: 



(1.21) 



Uz) a em (z) 



m 



where a lm denotes the components of the inverse of the matrix A = (a^ m ). 

The matrix A contains 0(h 2 ) terms that result from the pairing of directions paral- 
lel to the first soliton with directions parallel to the second soliton. Moreover, H(u z ) 



contains additional 0(h ) terms arising from the quadratic part of (1.14) not repre- 



sented in (1.19). It turns out that 0(h 2 ) terms in a im and 0(h 2 ) terms in H(u z ) each 



give rise to terms which cancel in (1.21). This hinges upon the fact that 



(1.22) 



where 



d y tH(u 7 



d Vl H(rn 



Vaje + d z e(vi, H' p {r) 2 )) + d z e(r) 2 , H' p (r]x)) 

3=1 



d V2 H(r] 2 



-d^H(rn) , b 7 = -d^Hfa) 



and all other bj = 0. When this equation is substituted into (1.21), once can witness 



the simplification arising from the pairing of A and A 1 , and this shows that (1.21) 



is equivalent to (1.19). We elaborate upon this in Appendix A. 



The merit in this point of view is that the equations (1.21) readily follow from the 



symplectic decomposition of the flow-that is, we select z (via the implicit function 
theorem) so that 



(1.23) u = u z + w 

where w G TzM 1 - (the symplectic orthogonal complement to T Z M in T Uz L 2 ) 



In 



(Lemma 3.1) a general argument is given showing that the equations (1.21) follow 
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with errors of size h A + ||w||#i 
errors of size h A + is equivalent to 



This argument exploits the fact that (1.21), with 



(1.24) 



d t u z = U z JH\u z ) + 0(h 4 + \\w\\ 2 Hl ] 



where U z : T Uz L 2 — > T Z M is the symplectic orthogonal projection operator given 
explicitly by 



n 2 /= ^(/.r 1 ^^)^™^)^ 



The proof of Lemma 3.1 makes no reference to the specific meaning of H or u z , and a 
similar result with nearly identical proof would yield the equations of motion in many 
other problems, including those studied in [HI [TJ [8] . The fact that the equations of 



motion follow automatically but rigorously from the symplectic decomposition (1.23) 



is one of the main advantages of this geometric approach to our problem, as opposed 
to a more ad hoc approach based on the discussion surrounding (1.16). 

It then remains to show that ||n;(i)||#i < h 2 on the time scale 0(/i _1 ), which we 
would like to prove using a suitable adaptation of the Lyapunov functional method 
initiated into the theory of orbital stability of single solitons by Weinstein [19] . Un- 



fortunately, the presence of the IT 2 projection in (1.24) corrupts this computation 



and only yields a bound ||u>(i)||#i < h. To eliminate this problem, we construct a 
function v z = 0{h 2 ), whose only time dependence is through the parameter z, such 
that the distorted double-soliton function u z = u z + v z satisfies 



(1-25) 



d t u z 



JH' 



0(h 4 + 



\w 



1 2 

Ih 1 



which is just (1.24) without the Il 2 projection. The construction of v z is carried out 
41 



in 



We add this correction v z to our soliton manifold M and consider the distorted 
manifold M = {u z }- The solution u to (1.1) now has a decomposition u = u z 



+ w 



where u z satisfies ( |1.25[ ) and it suffices to prove that < h 2 . In other words, 

we would like to show that the exact solution to (1.1) is approximately equal to 



the solution to the approximate equation (1.25). In £j5j a Lyapunov functional is 
employed to obtain the needed control on w. The Lyapunov functional used is a 
superposition of two copies-one for each soliton-of the classical functional, built from 
energy, momentum and mass, employed by Weinstein [19] to prove orbital stability 
of single solitons. This superposition was previously used by Martel-Merle-Tsai HI 



5 The idea that the equations of motions should be Hamilton's equations for the restricted Hamil- 
tonian with respect to the restricted symplectic form was introduced in [71 [8] and supported in- 
formally with an argument involving Darboux's theorem. The equations of motion thus obtained 
were used as a guide in the analysis in [51 [71 15] but the general rigorous connection between the 



symplectic decomposition of the flow and the equations of motion, as obtained in our Lemma 3.1 
was not obtained in [H [7J [S] . 
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in their study of the orbital stability of spreading multiple solitons. Our presentation 
of this component of the argument is a little different from [TO] or [TO] and more in 
line with the abstract orbital stability theory developed by Grillakis-Shatah-Strauss 
[U E] . Roughly, we prove that if W : L 2 — > R is a (densely defined) functional such 
that the derivative is of order 0(h) on M, and if we define L to be the quadratic 
part of W above M, then dfL is essentially the quadratic part of the Poisson bracket 
{H, W}(u) above M, which we show is of order 0(h 5 ). 

Let us note that h~ 5 losses occur in several estimates, which were not necessarily 
indicated in the above introduction, owing to the fact that in the attractive case \v\ 
can exceed h by a factor of log h~ x and a decreases below ao, as well as the presence of 
an x-multiplication factor in terms involving d v .u z in both the attractive and repulsive 
cases. We indicate the presence of such losses by writing, for example, h 4 ~ . These 
losses are more carefully quantified in the concluding summary of the proof in §|6j 

We emphasize that the methods in ^3] and §|5j although stated only for the problem 
at hand, are fairly general and widely applicable to problems in orbital stability of 
single [TOJ IU E] and multiple [TT] solitons and the dynamics of solitons in slowly 
varying potentials [SI El El El H3], weak rough potentials [U El [TO], and the interaction 
of two soliton tails, as considered here. The portion of the analysis most specific to the 
problem at hand appears in £j4j where the approximate solution is constructed. In this 
section, we consider the two components of the double-soliton separately and exploit 
the group structure of each individual soliton to pull-back to a nearly-stationary 
problem, which can be solved by operator inversion. This method was introduced by 
Holmer-Zworski [8] to produce an improvement of the result by Frohlich-Gustafson- 
Jonsson-Sigal [3] on the dynamics of single solitons in a slowly-varying potential, 
eliminating the uncontrollable errors in the ODEs appearing in [3]. 

Let us point out some related papers. Marzuola-Weinstein [12] consider the dynam- 
ics of symmetric and antisymmetric states in a double well-potential. Krieger-Martel- 
Raphael [S] construct two-soliton solutions with separating components asymptoti- 
cally as t — > +oo for the nonlinear Hartree equation, where the long-range effects of 
the nonlinearity complicate the analysis but also lead to nonnegligible perturbations 
of the asymptotic trajectories. Our analysis is similar in several ways to that of [9], 
although our priorties are different - we study the dynamics for a finite (but dynam- 
ically significant) time of an initial data that is close to a double-soliton, wheras they 
provide infinite time dynamics for an exact double-soliton solution. The problem of 
stability of nonintegrable NLS multiple solitons, with components that separate as 
t — > oo, has been considered by Perelman [14] , Rodnianski-Soffer-Schlag [16] . and 
Martel-Merle-Tsai [TO]. 

We now remark on where we rely upon the even/odd symmetry assumption on 
the solution. While the arguments in ^yielding (1.21) apply in general, in when 
constructing the solution u z to the approximate equation ( 1.25 ), we do make use of the 
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symmetry assumption, although we have sketched an argument (not included in this 
paper) showing how one can adapt the argument to the general case. The symmetry 
assumption also greatly simplifies the computations carried out in Appendix A which 
ultimately yield the ODEs (1.9), (1.10), (1.11) in Theorem |1.1| The integrals in the 
general case appear very complicated, and we are less confident that we could control 
the propagation of 0(/i 4 ) errors, as previously remarked. However, the one place 
where the symmetry assumption is used critically is in obtaining the upper bound 
on the Lyapunov function used in $5] to show the closeness of the true solution u 



and the solution u z of the approximate equation (1.25). Our guess is that to resolve 



this issue, one would need to restructure the Martel-Merle-Tsai Lyapunov function 
in a substantial way. The lower bound on the Lyapunov function, however, carries 
through in general. 

1.1. Acknowledgements. We thank Maciej Zworski and Galina Perelman for help- 
ful discussion related to this paper. J.H. was partially supported by NSF Grant 
DMS-0901582 and a Sloan Research Fellowship. 

2. Background on solitons, Hamiltonian structure, and Lyapunov 

function als 



The NLS equation (1.1) can be put into Hamiltonian form as follows. Take as 



the ambient symplectic manifold L 2 = L 2 (IR; C) with metric {vi,V2) u = Re J V\V2 for 
u G L 2 , v±,V2 € T U L 2 = L 2 . Let J = —i, viewed as an operator T U L 2 — > T U L 2 . 
The corresponding symplectic form is u) u (vi,V2) = (v±, J~ 1 V2) U (we henceforth drop 
the li-subscript). Define the (densely defined, with domain D = H 2 ) Hamiltonian 
H : L 2 -> 



as 



H(u) 



\U\ 



1 

4 j ' — ' — 4 

Using the metric (-, •) defined above, H'{u) G T*L 2 is identified with an element of 



T U L 2 . The free NLS equation (1.1) is 
(2.1) 



d t u = JH'iu) 



Solutions to (1.1 ) also satisfy conservation of mass M{u) and momentum P(u) } where 

1 r, ,n __, , i 



\u\ 



Pu) 



Im uu x 



Let d>(x) 



M(u) 
sechx and 

r](x, fi, a, 6, v) = e ie e^' lv(x - a) fi(f){fi{x - a)) 

Direct computation shows that M(r]) = fi and P(i]) = v. Consider the manifold of 
solitons 

M = { rj(-, fi,a,9,v)\ n > 0,9 eR,a eR,v eR} . 
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The tangent space at r\ = r](-, fi, a, 9, v) is 

T( P) a,e,v)M = span{ d^rj, d e rj, d a rj, d v i] } . 
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Note that JH'(rj) G T^^^M, and thus the flow associated to (1.1) will remain 
on M if it is initially on M. Specifically, direct computation shows 



(2.2) 



JH'(v) 



'^-(i V + ^fi 2 )d e r] + /i 1 vd a T). 



To gain a better understanding of (1.3) and (2.2), we can restrict uj to M to obtain 

i*u = d6 A dfi + dv A da , 
where % : M — > L 2 denotes the inclusion and restrict H to M to obtain 

1 1 o 1 q 



H(rj) 



-fx "V 



6 



and then note that (1.3) is just the solution to the Hamilton equations of motion for 
H{rj) with respect to i*u: 



(2.3) 



( fi = d 6 H( V ) = 

a = d v H(rj) = fjT l v 
1 



-22,1 2 



6 = -d^H( V ) 
, v = -d a H( V ) = 

Suppose we knew that JH'(rj) G T^^q^M and wanted to recover the coefficients as 
in (2.2). This could be achieved by noting that 

JH'(rj) = d t rj 

= fid^r] + ad a 7] + 8d e rj + vd v rj 
= d v H(rj)d a r] - d,,H(rj)derj 

Moreover, the functionals M and P, considered as auxiliary Hamiltonians, have as- 
sociated Hamilton vector fields 

JM'(r ] ) = -d e r ] JP'(r ] ) = d aV . 

This enables us to write 

(2.4) JH'(rj) = d v H(rj) JP'(rj) + M #(ij) JM'fa) . 

0, where 



From this, we learn that W'^ a & v {rj) 



(2.5) 



W, 



(u) = -d^H(rj)M(u) - d v H{ v )P{u) + Hi 



U 



The functional L^ afi ^ v) (u) = W^ a fi, v ){u) 
used in the classical orbital stability theory for (1.1) due to Weinstein |19j . 



W(/x,a,0,*O (^) ^ S ^ ne Lyapunov functional 
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3. Effective dynamics 



Now we turn to the double soliton problem and begin the proof of Theorem 1.1 
Consider the two-soliton submanifold M of L 2 given by 

M = {u z = r](-,fi 1 ,a 1 ,6 1 ,v 1 ) + r)(-,/i 2 , a 2 , 9 2 ,v 2 ) } . 

Note that M is just the linear superposition of two single solitons. We adopt the 
notation 

z = 0\ z 2 , z 3 , z 4 , z 5 , z 6 , z 7 , z 8 ) = (/ii, ai, \i 2y a 2 , 9i, Vi,0 2 , v 2 ) , 

for coordinates on this manifold M. Next, we give the form of the symplectic orthog- 
onal projection operator 

U Z :T U L 2 ^T Z M, 

Note that T Uz L 2 is naturally identified with L 2 . A consequence of the requirement 
that (/ - ILJ, J~ x d z iu z ) = 0, I = 1, . . . , 8 is that 

8 

(3.1) U z f = J2 (/> J~ 1 d z m z }a em (z)d zmUz 

£,m=l 

where A(z) = (a£ m (z)) is the 8x8 matrix with components ae m (z) = (d z m z , J~ x d z ™u z ) 
and A(z)~ l = (a em (z)) is the inverse matrix. 

Let i : M — » L 2 denote the inclusion. It follows from the definition of A(z) that 
the restricted symplectic form takes the form 

1 8 

(3.2) i*(u) = - a ^dz l A dz m . 

£,m=l 



It also follows by substitution into (3.1) that 

8 

U z JH'(u z ) = -J2 d z iH{u z )a tm {z)d zmUz 
e,m=i 

Consequently, the equation d t u z = U z JH'(u z ) is equivalent to the system of equations 

8 

z m = - 9 z eH(u z ) a £m (z) m = l,...,8, 

which are precisely Hamiltonian's equations of motion for the restricted (to M) Hamil- 
tonian z h- > H(u z ) with respect to the restricted (to M) symplectic form i*(u>). 



We propose to model the solution u to (2.1 ) by 

(3.3) u = u z + w 

where u z G M is chosen so that the symplectic orthogonality conditions 

(3.4) (w,J' 1 d z m z )=0, £ = !,..., 8. 
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hold. The fact that such a z exists follows from the implicit function theorem and 
the assumed smallness of w. Note that if we assume u(t) solves (2.1), this induces 
time dependence on the parameters z £ K 8 . 



Lemma 3.1 (effective dynamics). Suppose that u evolves according to (2.1) and z, 
w are defined by ( |3.3 ) so that the orthogonality conditions (3.4) hold. Then 

(3.5) \\d t u z - U z JH'(u z )\\ TzM < \\w\\ 2 m + max || J- l d z ^JH'{u z )f m . 

l<n<8 

Equivalently, considering M as an 8- dimensional symplectic manifold equipped with 
the symplectic form i*(u) given in (3.2), the Hamilton's equations of motion for z 
induced by the restricted Hamiltonian z i-> H{u z ) approximately hold as follows: 
(3.6) 



z m 



+ J2^H(u z )a em (z) 



< \\w\\ 2 m +msx \\J^d zn UUH'(u z )\\ 2 m 



Kn<8 



m 



The norm || • \\t z m is the one induced by the metric (•,•)«*■ As T Z M is finite- 
dimensional, we have the norm-equivalence to 



^2j(z e )d z m, 



~ max \^(z )\ 
i<e<8 



Proof. Since u solves (2.1), we obtain from (3.3) the equation for w. 

(3.7) d t w = -(d t u z - U z JH'(u z )) + UjjH'(u z ) + JH"(u z )w + O m (\\w\\ 2 Hl ) 



By applying dt to (3.4), we obtain 

= (d t w, J' x d z -au z ) + (w, J' x d z r,d t u z ) 

Here we have used that <9f<9 z n = d z ndt, which holds provided we adopt the con- 
vention that d z ez m = for all 1 < £,m < 8. Substituting (3.7) and using that 
(n. z JH'(u z ), J~ 1 d z ™u z ) = 0, we obtain 



(3.8) 
where 



0=A+B+C+D 

A = -{d t u z - U z JH'(u z ), .r x d z ^u z ) 
B = (JH"(u z )w, r l d zn u z ) 
C = (w, J~ l d z nd t u z ) 
D = (0(w 2 ), J~ x d z nu z ) 



Note that here w is properly understood as an element of T U ^L and in (3.3) we mean that, 
starting at u z we take the flow- forward (by "time" 1) in the direction w. However, using the natural 
identification between T Uz L 2 and L 2 , (3.3) makes sense as an equation involving functions in I?. 
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Since J* J~ l = — 1 and H"{u z ) is self-adjoint, 

B = —(w,H"(u z )d z nU z ) = —(w,d z nH'(u z )} = —(w,J~ 1 d z nJH'(u z )) 

Hence 

B + C = (w, r x d zn [d t u z - JH'{u z ))) 

= (w, r l dAd t u z - u z jh'(u z ))) - (w, r l d z ^jH\u z )) 

Let R = d t u z — H z JH'(u z ) G T Z M, and expand with respect to the basis of T Z M as 



R = ^2j e (z)d z 



It follows from (3.8) that 

(3.9) (R, J~ x d z -au z ) = (w, J~ l d zn R) - (w, J^d.nUj JH'{u z )) + 0(\\w\\ 2 Hl ) 
We have 

8 8 

d zn R = ^ 9 z nji(z) d z tu z + ^ le(z)d z nd z m z . 



Since w G T Z M^, 



J 1 d z nR) = ^2>y e (z){w,J 1 d z nd z m z ) 



i=i 



and hence 
(3.10) 



{w^-^nR}] < \\w\\hi\\R\\t, 



M ■ 



The lemma follows from (3.9), ||-R||t z m 
Cauchy-Schwarz. 

In our case we shall have 



maxi<„< 8 | (R, J 1 d z nu z )\, (3.10), and 



□ 



\\J- l d z ^JH'(u z )\\ 2 m <h 4 - . 

We carry out computations of (3.5) in Appendix |A| and show that (3.5) is equivalent 
to (1.19), with error terms 0(h 4 ~), even without the even/odd assumption on the 
solution. It is further shown in Appendix |A] that when the even/odd assumption is 
imposed and the integrals in (1.19) are explicitly computed, we obtain 

' (i = (-l) CT (8a - A)ve~ 2a + 0(h 4 ~) 
a = yT x v + (-l) CT (-4a + \ 2 )ve- 2a + 0(h 4 ~) 

0=l^ 2 + ^V 2 + 18(-l) CT e- 2a + 0(h 4 ~) 

v = -4(-iy e - 2a + 0{h 4 -) 
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The solution (p, a, 9, v) is adequately approximated by the ODEs appearing in the 
statement of Theorem 11.11 

4. Approximate solution 



By Lemma 3.1 and (3.7), the equation for w is 

= JH"{u z )w + UjjH'(u z ) + O h x{ 




\w\ 



H 



x + h 4 



t=o 



w 



The next step is to show that there exists a function v z (x) such that < h 2 , 

whose only time dependence occurs through the parameter z, such that 

(4.1) d t v z = JH"(u z )u z + TljjH'(u z ) + H x(h 4 - + \\w\\ 2 m ) 
Here it is assumed that zeK 8 evolves according to Lemma 3.1 , i.e. 

(4.2) d t u z = Tl z JH'(u z ) + H i(h A ~ + \\w\\ 2 H1 ) 

The initial data ^L_ is not prescribed but our structural assumption on v z is fairly 
rigid. Note that given (4.2), the assertion that v z solve (4.1) is equivalent to the 
statement that u z == u z + v z solve 



(4.3) 



d t u z = JH\u z ) + O m {h 4 ~ + \\wf m ) . 



This is an approximate solution to (2.1 ) (that does not, in general, satisfy the specified 
initial data). 

Let g : L 2 — > L 2 be the operator attached to the parameters (fi,a,8,v) that acts 
on a function p as follows: 

(4.4) (gp)(x) = e w e^^pp^x - a)) . 

The inverse action is 



g 1 p{x) 



-Ifl "vx 



p 1 p(/i x x 



The adjoint action g* with respect to (■, ■) is 



g*p(x) 



-Ifl "vx 



p(p x x + a) = pg l p(x) . 



Denote 4>(x) = sechx. Then u z = g\<\> + g2<p- We look for a solution v z to (4.1 ) in the 
form 

2 

(4.5) u z = J2a jgjPj , 

i=i 

where <x,- = acj(pi, a%, 8i, Vi, P2, ^2) ^2) an d gj is the operator corresponding to 
(pj,aj,9j,Vj). That is, we assume v z can be decomposed into two pieces, each of 
which can be pulled back to a stationary equation and solved by operator inversion. 
The time dependence of v z occurs only through z. 
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The next step is to substitute (4.5) into (4.1). The resulting equation simplifies 
provided we assume that each pj satisfies < h 2 e^~ 1+ '^ for \x\ > 1 as certain 

cross terms become 0^i(/i 4 ~). In this case, (4.1) will be satisfied provided for both 
j = 1, 2, we have 



dtiajgjPj) = JH"(g j (p)(a j g j p j )+U z J((g j (f)) g 3 _ j <f)+2\g j <j)\ g 3 ^(j))+0 H i(h +\\w\\ H i) 

In the proof, we delete the j-subscripts, denote g = g 3 ^j and moreover assume that 
ocj = 0(h 2 ). Then we aim to solve 

d t (gp) = JH"(g^)(gp) + a^U 1 J^fM + 2|#| 2 #) + O m {h A ~ 



\w\ 



The form of the operator IT- 1 - can be simplified, since we only need to keep the 0(1) 
and 0(h) parts. This equation takes the form 

(4.6) d t {gp) = JH"(g^)(gp) + a'^Jgf + O m (h 4 - + \\w\\ 2 m ) 

where, in the case j = 1, 



(4.7) 



/ =srMM) 2 ^] + %r%i0iW] 

= ^^p 2 p 3 ^{p 3 ^pfx + ( aj - a^))<P(xf 
+ 2e^e^ x p 2 p^ j( f>(p 3 ^pJ l x + ( aj - a 3 ^))4>{xf 



with 



dcf 



0i - 9 2 -/i 3 i f 3 -i(aj 



a 3 _ 



dcf _2 -1 -1 

W 2 = Pj Vj - Ps-jPj j 



U) 3 



dcf 



In the case j = 2, 
(4.8) 



dcf _2 -1 -1 

W 4 = Pj Vj + P^jPj j 



with a similar expansion. The only important feature of these expressions is that 
/ = 0(h 2 ) and e^~'^dtf = 0(h 3 ) when 9\ — 2 is a constant (as in the case 



,(!-)<*) 



of the even/odd symmetry assumption in Theorem |1.1| 

Now we begin the task of pulling back (4.6) - applying g* to (4.6), we obtain 

(4.9) g*d t gp = g*[JH"{g4>){gp)] + a^g^Jgf + O m (h 4 - + \\w\\ 2 m ) . 



First, we aim to simplify the term g*[JH"(g<j))(gp)] in (4.9). Let K g (4>) = H(gcj)). It 
follows that K' g (<j>) = g*H'(g<j)) and K'^) = g*[H" (g4>)(g»)\. By direct substitution, 
we compute: 

K g (u) = p 3 H(u) + pvP{u) + \v 2 p~ l M(u) 
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Since g* J — Jg*, we have 

g*[JH"{g$){g.)] = JK>>{4>) 



(4.10) 



H 3 JH"{<f)) + (ivJP"{<f)) + \v 2 iT l JM"{(j)) 



Second, we seek to simplify the term g*d t gp in (4.9). Define the operators 





def 




(1,0,0,0) 


d a 


dcf 


da 


(1,0,0,0) 


Oe 


def 


de 


(1,0,0,0) 


d v 


def 


d v 


(1,0,0,0) 



d x x 

—% 

i 

ix 



Let 



(4.11) 



def * 

o» = g o^g 

}\ def * r% 

o a = g d a g 

k def * 

de = g o e g 

p. def * 

d„ = g d v g 



d x x — in 2 vx 
-p 2 d x - iv 

iji~ l x 



<9 M - p~ 2 vd v 
li 2 d a - vd e 
lid e 
H~ l d~v 



It follows from the chain rule that 

g*d t g = jjBp + ad a + 6d e + vd v . 



Using (3.5) 
(4.12) 



g*d t g = vpL 1 d a + \^d e + H i(h 2 
= v/jd a + y 3 d e + H i(h 2 ) 



Finally, we aim to simplify the term g*U ± Jgf in ( |4.9 ). We will show that 

(4-13) 9*K(ji,a,e,v)Jgf = ^11(1,0,0,0)^/ 

Using that J*J _1 = —1 and (4.11), we obtain 

9*ttJgf = (Jgf, J~ 1 d a g<p)g*d v g<f) - (Jgf, J- 1 d v g<f))g*d a g<j) 

+ (Jgf, J~ 1 d fl g<f))g*d e g<f) - (Jgf, r 1 d e g$)g*d^ 
= ~(f, d a 4>)d v (t) + (/, d v <p)d a <f) - (f, B^dgcj) + (f, de&d^ 



Substituting (4.11), after a few cancelations we obtain 

yT^UJgf = -(fA4>W+(fMW-(fA^W+(fM)d^ 

= + (Jf, J~ l da<P)d v <P - (Jf, J- l d v <P)da(t> 

+ (Jf, J~%4>)de<l> - (Jf, J^dgfldrf 



which establishes (4.13). 
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Note that it follows from (4.13) that 



9*R(n,a,6,v)9Jf — (^(1,0,0,0) J f 



(1,0,0,0)' 



Using the expressions (4.10), (4.12), and (4.13), the equation (4.9) converts to 



vpd a p + \p z d e p + pd t p = fi 3 JH"((P)(p) + fivJP"(4>)(p) 

+ a-Vnf li0 , ,o)^/ + O m (h 4 - + || 
Noting that JP"(4>) = d a and JM"(0) = — dg, the equation becomes 
p 3 JM"((f>)p + p?JH"{<j))p - pd t p = -a-Vnf Jf + O m (h 4 - 



H 1 



1 ,.3 

2 



\W\ 



Hence we see we should take a = p 2 so that the equation becomes 

\JM\<\>) P + jh'\<\>)p - p- 2 d tP = -nf li0)0)0) Jf + o m (h A - + \\w\\ 2 m ) 

Now apply J -1 to obtain the equation 

(4.14) S P = r l p- 2 d t p - r l Tlf lflm Jf + 0{h 4 ~ + \\wf m ) 
where the operator 

S(p) d = f §M"(0)(p) + tf"(0)(p) = \p- \d 2 x p - 2|0| 2 p - 2 p 

is self-adjoint with respect to the inner product (u, v) = Re J uv. The kernel is 
spanned by d a <t> and dg(j). 

Lemma 4.1 (properties of S). 

(1) For any f G H 1 , let F = J -1 !!^ ^ Jf. Then F satisfies the orthogonality 
conditions 

(4.15) {F,d e <t> 

(4.16) (F,d^ 



0. 

0, 








(2) For any F satisfying (4.15), S l F is defined and satisfies the boundedness 
properties 

(4-17) WS^FWhx < \\F\\ L 2, 

(4.18) We^S^FWnz < ||e ff(:c> F|| L 2 . 

for < cr < 1 . 



(3) For any F satisfying (4.15) and (4.16), S l F satisfies the orthogonality prop- 
erties 



(4.19) 
(4.20) 



(Sr L F,d 9 <l>) ~- 
(J-'S-'F^e 



0, 



(S- l F,d a 

(J-'S^F^a 
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Proof. Item (1) is immediate from the definition of 11(10,0,0)- F° r item (2), we recall 
that ker S = span{<9 a 0, de<f\ and moreover, is an isolated point in the spectrum of 
S. Thus S^ 1 : (kerS 1 ) -1 — > (ker S) 1 - is bounded as an operator on L 2 . The inequality 
(4.17) follows from this and elliptic regularity. To prove (4.18), it suffices to show 
that for any G and any \a\ < 1, we have 

(4.21) \\G\\ H , < \\e™Se-™G\\ L 2 + He^GU^ 



Indeed, (4.18) follows by taking G = e ax S F, appealing to (4.17), and separately 
considering a > and a < with |a| < 1. To establish (4.21), we calculate 

(4.22) e ax Se~ ax G = (S + ad x - \a 2 )G 



and hence 



_ J) + adx _ \ d l )G _ 2 <p 2 G - <fG 



[I - a 2 ) + ad x - -d 2 x )G = e ux Se- ax G + 20 2 G + 2 G 



v 2 v * 2 

On the left-hand side, we have an operator with symbol |(1 — a 2 ) + ai^ + £ 2 , which 
dominates (£) 2 under our assumption on a. From this and the fact that 



[x 



< 



-2\x\ 



we conclude (4.21). 



For item (3), (4.19) follows from the fact that 5* 1 : (kerS 1 ) -1 — > (ker,!?)- 

(S- l F,Sd^) 



To 



establish ( |4.20[ ), we note that by ( |4.16[ ) 

0=(F,d^ 



and similarly 

= (F, d v 4>) = (S^F, Sd v <f>) 

and thus it suffices to establish that S(d v (p) = J d a $ and S'((9 M 0) = J~ 1 dec 
prove these equalities, recall 

(4.23) 

Taking d v and evaluating at (/x, a, 9, v) = (1, 0, 0, 0) gives 

S(d v <f ) ) = P'(<j ) ) = J~ 1 d a <p. 
Taking of (4.23) and evaluating at (fi,a,9,v) = (1,0,0,0) gives 

Sid^) = -M\<t>) = J- l d ef i 



To 



= W'(rj) = (i/i 2 + V)M'(7?) - ^vP'irj) + H'( V ) 



□ 



Recall that the task is to solve (4.14) where / is either ( |4.7[ ) or ( |4.8[ ). At this point, 
we impose the even/odd solution assumption as in Theorem 1 1 . 1 1 which implies that 
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Q\ — 8 2 is constant. The other time dependent parameters in (4.7), (4.8) are all slowly 
varying so that d t fj = 0(h 3 ). Thus, we can solve (4.14) by iteration. LetQ 



(4.24) 



9 



—S l j 1 n( 100)0 ) Jf 



By Lemma 4.l[ l)(2), this is well-defined with p 1 = 0(h 2 ~) and satisfying all the 
needed regularity properties. With p 2 as yet undefined, we plug p 1 + p 2 into (4.14) 
to obtain 

s P 2 = r x p- 2 d tP x + j- V 2 <V + o(h 4 - + HIIO 

As mentioned previously, d t f = 0(h 3 ~) and thus 



is also 0(h 3 ). By Lemma 4.1^3), in particular (4.20), with F = J 1 I1^ L 1 ^ Jd t f, 
we have that 

F = J^p-^p 1 = -J- x p- 2 S- l F 

satisfies the condition (4.15), and hence we can apply Lemma 4.1| (2) with F replaced 
by F. That is, the function 

(4.25) p 2 d ^ S-'F = -S-Kr^-'S-'J-'Uf^JdJ 

satisfies all the needed regularity properties. Note further that d t p 2 = 0{h A ~). Upon 
substituting p 1 + p 2 into (4.14) with p l defined by (4.24) and p 2 defined by (4.25), we 
find that equality holds with 0(h 4 ~) error. 

Thus we have successfully constructed a solution to the approximate equation (4.1 ). 
We summarize our conclusions in the next lemma. 



Lemma 4.2 (approximate solution). Recall the operator gj associated to (pj, aj, 6j, Vj) 
defined in (4.4) and fj defined in (4.7) (j — 1) or (4.8) (j = 2). Let p l - be given by 
(4.24) and then let p 2 be given by (4.25). Then p k - for 1 < j, k < 2 satisfy 



\e {1 - ){x) p k Am < h 



j 

l+k- 



Let 



Vz = A_9\{.9\ + 9 2 i) + M2^2(P2 + p\) 

Suppose that the parameter z G M 8 evolves according to the ODEs obtained from 
Lemma 3.1 (in the same phase or opposite phase case). Then v z {x) solves (4.1). 



^As indicated earlier, p can stand for either p\ or pi. The superscript introduced here is different 
and meant to indicate part of an asymptotic expansion for either function. In other words, we have 
Pj = p) + p] for both 3 = 1, 2. 
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5. LYAPUNOV FUNCTIONAL 



The final step is to show that the true solution u to (2.1) is approximately the 
approximate solution u z = u z + v z . For this purpose, we introduce a Lyapunov 
functional. First, some general considerations. We consider the "perturbed" 8- 
dimensional manifold 

M = { u z | z G M 8 } 



Introduce the notation w = u — u z (so that w = w + v z ). Now it follows from (4.1) 
that 

(5.1) d t u z = JH'(u z ) + F. 

where F = H i(h A ~ + ||u>||^i). 

Suppose that W z : L 2 — > E is a densely defined functional. We write d z eW z : L 2 — > 
M. to indicate partial derivatives with respect to z and 

W' z {u) G T*L 2 ~ T U L 2 ~ L 2 

metric g 

to indicate partial derivatives with respect to it (ignoring the interdependence between 



z and u given by (3.3), (3.4)). 

Suppose that W z can be extended to a differentiable functional H 1 — > K; then 

for each n G i? 1 , we have a bounded linear map W' z {u) : ff 1 — > R which, under 

the aforementioned identification, becomes a function belonging to H~ x . In fact, our 

choice of W z is differentiable at all orders as a map H 1 — > IR, which is to say that 

x ■ ■ ■ x H 1 — R is a bounded /c-multilinear map. 
s v ' 

fc copies 

We further assume that d z iW z (u) = unless \z \ < h. Let 
(5.2) L z (u) = W z (u) - W z (u z ) - (W& g ),w) . 

That is, L z (u) is the quadratic part of W z (u) above the base manifold M. Now 



viewing u = u(t) and z = z(t) in accordance with (3.3), (3.4) (and thus reinstating 
the interdependence between z and u), we have, for any functional G z : L 2 — > M, 

8 

d t G z (u) = (G' z (u),d t u) + J2 [9** G] {u)z k . 

k=l 

This leads to: 



Lemma 5.1. Suppose that u solves (2.1) and z evolves so that u z solves (5.1), and 



that L z {u) is given by ( 5.2[ ), the quadratic part ofW z (u) above M. Then 



(5.3) d t L z (u) = {H, W z }(u) - {H, W z }(u z ) - ({H, W z }'(u z ),w) -E l + E 2 
where 

El ^ (W"(u z )F,w) + (W'(u z ), \JH\u) - JH\u z ) - JH"(u z )w}) 
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and 

8 

E 2 * J2([d^W z ](u) - [d zk W z ]{u z ) - ([d zk W' z \{u z ),w))z k 

k=l 

In other words, d t L z (u) is, up to error Ei and E 2 , the quadratic part of {H,W z }(u) 
above M. Note that E 2 just involves the quadratic part of [d z kW](u) above M. 

In the typical application of this lemma (as for our W z , defined below), we have 
bounded operators W z {u z ) : H 1 — > H^ 1 and W z "(u z ) : H 1 x H 1 — >■ L 2 which implies 
the bound 



\Ei\ < \\F\\ H i\\w\\ H i + \\W'(u z )\\ H i\\w\\ 2 H i 



Thus, one just needs < h 3 and || W'{u z ) \\ H i < h; in our case we in fact have 

the stronger statements < h 4 ~ and \\W {u z )\\ H i < h 2 ~ . Moreover, in our case 

we will have 



<h 2 ~\\w\\ 2 Hl 



since fi,v = 0(h 2 ). 



Proof. By (5.2) 



d t L z {u) = d t W z {u) - d t W z {u z ) - d t (W' z {u z ),w) 
We compute each of the three terms on the right-hand side separately. 



d t W z (u) = (W' z (u),d t u) + Y^[d^W z ]( 

k=l 

8 

(5.4) = (W' z (u), JH\u)) + J2l^W z ](u)z k 

k=l 



where we invoked (2.1). Second, we compute 

8 

d t W z (u z ) = (W' z (u z ),d t u z ) + Y} d * kW *\( 

k=l 

(5.5) = (W z (u z ),JH\%)) + (W z (u z ),F) +Y^[d z «W z }( 



• )z k 

/,-: J 

8 

k 



.^U^Z 

k=l 



NLS SOLITON INTERACTION 



2:5 



where we invoked (5.1). Finally, we compute 

d t (W' z (u z ),u-u z ) 

= (W"(u z )d t u z , u-u z ) + (W' z {u z ), d t u - d t u z ) 

8 

+ Y,([dz*W}'(u z ),u-u z }z k 

k=l 

(5.6) = (WJ{u z )JH'(u z ),w) + (W"(u z )F,w) + (W' z (u z ),JH\u) - JH\u z )) 

8 

- (W'(u z ),F) + (J2[d zk W]'(u z ),w)i k 

k=l 

Taking Q minus ( |5.5[ ) minus ( |5.6[ ), noting the cancelation of +(W' z (u z ), F) in ( |5.5[ ) 
with -(W'(u*), in Q, we obtain □ 

To produce W z (u), we use an idea of Martel-Merle-Tsai [ID] . Let 

1 if x > 1 



if x < — 1 



such that (^ / ( a; )) 2 ^ min(*(a;), 1 - *(^))- Set $ = 4 /( lo g = 4 / a o, so < 5 < 1. 
Introduce the localizations ^2 (2) = ^(foc) and t/>i(:e) = 1 — ty(5x), and set Mj(u) = 
M{i)) /2 u) and Py(u) = P{^) /2 u). Define 

(5.7) d = f - ± d -^M A u) - ± ^PAu) + H(u) 



^2 2 



i=i 



The Lyapunov functional L z (u) we use is then defined as in (5.2). 

Lemma 5.1 facilitates the computation of dtL z (u), since W z (w) is built from "nearly 
conserved" quantities. Indeed, we have the following Poisson brackets: 

{H, Mj}(u) = \^J $'juu x 
{H,Pj}(u) = J V^W 2 W 



'■"'\u\ 2 



It thus follows from Lemma [5. II that 
(5.8) 



d t L z {u) =-^(//f + ^\f)({#,M i } // (^ 



3=1 



+ 0(\\w\\ 3 m )-E 1 + E 2 
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For our choice of W z (u), as remarked earlier, we have suitable bounds for E\ and E 2 . 
Moreover, once one imposes the even/odd solution assumption of Theorem we 
have \i\ = fi2 and v\ = v 2 , so the first term in ( 5.8[ ) disappears]^] Hence 

(5.9) \d t L z (u)\ < ((K| + \v 2 \)5 + h)\\w\\ 2 m + h 3 \\w\\ Hl 



7.113 



W 



\H 1 



Since \vj\ < h 1 \ogh 1 and 5 ~ (log/i x ) 1 , the term (\v\ \ + |t>2|)<5 < ft<- 

Now we turn to the matter of obtaining a lower bound for L z (u). First note that 

(W' z '(u z )w,w) = L z (u) + 0(\\w\\ 3 Hl ) . 

Given that II u z \\ h 1 ^ h 2 ~ , we have 



(5.10) 



(W'J(u z )w, w) = L z (u) + 0(h 4 -) + O^IMI^ . 



The needed lower bound for the left-hand side will be established below in Lemma 

E2 

For the single-soliton case, we have coercivity for the classical functional from 
Weinstein [TH], which we now recall. Taking rj = fi, a, 9, v) and 



R 



def 



dH{rj) 
djj, 



M{u) 



dHjjj) 
dv 



P{u) + H(u) 



(/i 2 + /j- 2 v 2 )M(u) - jj^vPiu) + H(u) 



then 
(5.11) 



< 



(R"(r])w } w} } 



provided we assume the orthogonality conditions 
(w } J- 1 d flV } = 0, 
(w, J-'dev) = , 



(w, J- l d aV ) = 0, 

(W, J^dvT]) = . 



A direct proof of (5.11) is possible; see [7, Prop. 4.1]. 

We now prove a similar argument for the double-soliton functional W z (u) defined 
in (5.7). Before proceeding, we record the formulae 

M"(u) 



(5.12) 



H"(u) 



l -d 2 

2 x 



2\u\ 



u 2 C 



where C denotes the operator of complex conjugation. 



In fact, this is more easily seen by observing that once [i\ = \ii and v\ = v%, we have that the 
first term in (5.7) becomes M(u), whose Poisson bracket vanishes. We included the localization in 
this term to illustrate the difficulty in treating the asymmetric case - one would not have that the 
first term in (Ol is 0(h 5 ). 
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Lemma 5.2. Suppose w satisfies the orthogonality conditions (3.4). Then 

\\Hm<(W';{u z )w )W ). 



(5.13) 



1/2 

Proof. Denote Wj = r/y w, j = 1, 2. Note that u>i + w 2 ^ w, although 1 = + i\) 2 < 
ipy 2 + ^2^ 2 < 2. Define functional 



-Af(u) 



P(u) + if(u) 



r(/i 2 + fxj 2 v 2 )M(u) - fxfvjPiu) + H{u) 



We claim that 
(5.14) 

and 
(5.15) 



(W"(u z )w,w) - ^W^iH-^i) 



<5 2 llw; 112 



hi 



hi 



E 

3=1 



<5 2 |b" 2 



u 1 



We now establish ( |5.14| ). Note that (see Q5.12D ) 

2 



(w"(u z )w,w) - w j> 



3=1 

= ((#>,) - ^H"^)^ 3 - 4 /2 H"( V2 )4 /2 )w,w) 

The operator appearing on the right-hand side can be decomposed into Ai + A 2 + A 3 
where 

A 2 ^-2(H 2 -ViM 2 -^l 2 ) 

We compute Ai explicitly: 

2 2 
3=1 3=1 

where we have used that ^1 + ^2 = 1 m the second equality. We have (^) 2 < f^ 2 ^- by 
the corresponding property of \I/ and thus Ax is a multiplication operator with symbol 
bounded by 5 2 . By the support properties of ipi, ip2, we obtain that the multiplication 
operators A 2 , A3 have symbols bounded by h. This completes the proof of (5.14), 
and the proof of (5.15) is similar. 
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By the orthogonality conditions (3.4), 



(w u J l d zJ u z ) = -((1 - ip\ /2 )w, J l d z jU z ) , 
and we have, for example 

( Wl ,j-% lUz ) = -((l-^Kr^u,) = -((l-^V^ 1 ^) £ h 1/2 \M\v , 

due to the fact that ||(1 — ip^ 2 ) J^d^rji < h 1 / 2 . Hence, by the coercivity of the 



classical Lyapunov functional (see discussion surrounding (5.11)), we have that 



\w 



jUH 1 > 



i=l 



i=i 



From this and (5.14), (5.15), we obtain (5.13). 



□ 



6. Conclusion of proof 



In this section, we conclude the proof of Theorem 1.1 



Recall that h = e a ° which implies that = log h 1 , and that we are in the 



even/odd solution setting with (1.7), (1.8) in place. 

We introduce < (S < 1. The constant 5 is absolute and is chosen sufficiently 
small in terms of the accumulation of numerous other absolute constants appearing 
in several estimates. In our argument, c will represent a large absolute constant that 
may change (typically enlarge) from one line to the next. At the conclusion of the 
argument, we can finally declare that 5 should be taken small enough that c5 < |. 
This does not constitute circular reasoning since one could tally up all of the absolute 
constants (the c's) in each estimate in advance of executing the argument and suitably 
define 5 a priori but this is not a practical manner of exposition. 

Recall that we started by defining 

w(t) = u(t) - u z{t) 

where z was selected by the implicit function theorem so that orthogonality conditions 



(3.4) hold. By continuity of the flow in H 1 , this is possible at least up to some small 
positive time. Let T be the supremum of all times < T < h^ 1 ^ 5 for which 

|u>Hl°° m < h 3 / 2 
v\ < h 1 - 6 

1-6 



(6.1) 
(6.2) 
(6.3) 



a > a 



o 



Note that the requirement (6.1) implies 



Mini - h \\ w \\m + h 3 \\w\\ H i 



and enables us to discard cubic error terms in w in our estimates. 
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In the course of the argument that follows, we work on the time interval [0,T]. At 
the conclusion of the argument, we are able to assert that either T = 5h~ x \ogh~ 1 or 



that (6.2) or (6.3) fail to hold at t = T. 



It follows from the decomposition (1.22) and the bootstrap assumptions (6.2), (6.3) 
above that (see Appendix [A| 



Let 



sup WJ^d^JH'iu^WnxKh 

Kn<8 



def , A—rr) n n 2 

e = h 4 cd + y oo „i 

[O.T]-"^ 



2-c<5 



By Lemma 3.1 and the computations in Appendix [A] the ODEs 

' ft= (_i)-(8a-4)t;e- 2a + 0(e) 

a = ^~ l v + (-l)°(-4a + \a 2 )ve- 2a + 0(e) 

= l^ 2 + ^V 2 + 18(-l) CT e- 2a + 0(e) 
v = -4(-l) CT e - 2a + 0(e) 



hold on [0,T\. By the first of these equations and (6.1), (6.2), (6.3), we have \fi — 1| < 
h 2 ~ cS . From the above ODEs and JO, ([6~2|, (B, we can deduce bounds on fi, a, 9, 



and v that justify the estimates involved in the construction of v z in §4 summarized 
in Lemma 4.2. The result is that 



(6.4) 



\ v z\\m 



< h 2 ~ cS 



and (5.1) holds with 



By (5.9) 



(6.5) 



< e. 



\d t L z (u)\<h\\w\\ 2 m +h 3 \\w\\ m + 



~ 1 1 3 



'W 



where we recall that w = w — v z . By (6.4), < II^IIh 1 + h 2 cS , we obtain from 

( |6~5| ) that 

(6.6) \d t L z (u)\ < h\\w\\ 2 Hl + h 5 ~ c5 . 



From (5.13) and (5.10), the bound 

\\w\\ 2 m <L z (u) + h A - cS 



(6.7) 



holds. Combining (6.6) and (6.7), we obtain the bound 

\d t L z {u)\<hL 2 (u) + h 5 ~ cS 
By Gronwall's inequality, it follows that 



LJu) < e cm h 



cth >4— c<5 
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Provided we restrict to t < Sh^ 1 logftr 1 , this implies 

LJu) < h 4 ~ cS 



Reapplying (6.7), we obtain 
(6.8) 



,2-c<5 



IIHItf 1 < h 

At this point, we can declare that 5 should have been taken sufficiently small so that 



c5 < |, where c is as it appears in (6.8). It follows that (6.1), (6.2), (6.3) can only 



break down provided T > Sh 1 log ft 1 or if either (6.2) or (6.3) fails at t = T 



We will see the from the following ODE analysis that (6.3 ) always holds; in the same 
phase (even solution, attractive) case, the assumption (6.2) first fails at T ~ ft, -1 , and 



in the opposite phase (odd solution, repulsive) case, (6.2) remains valid and we can 
reach T ~ ft -1 log ft -1 . 



Since we now restrict to t < 5h 1 log ft x , we can assume that (6.8) holds and thus 
e<ft 4 " c5 - 

Let z = (/i, a, 9, v) solve 

{-iy{%~a-A)ve- 2 ~ a 



f ~ 



a = v 
1 



1 



= -jl 2 + -v 2 j2~ 2 + 18(-l) CT e -25 
'v = -4(-l)V 



-2(7 



These tilde equations appear in the statement of Theorem |1.1| without tildes. Note 
that the a and v equations can be solved separately as discussed in ^ Let a = fia — a 
and v = v — v. Then we get the system 

= v + 0(h 3 ~ c5 ) 
= -2(-l) CT e -2(i a + 0(h 4 ~ c5 ) 

Let 7 = (a) 2 + h~ 2 v 2 . Then, substituting 

7 < ha{hT l v) + (ft 1/2 a)ftl -c<5 + (h~ 1/2 v)(hl- cS ) 

By the inequality a/3 < a 2 + /3 2 , we obtain 

7 < ^ + h 5 - cS 

By Gronwall's inequality, 

7 < e cM (7o + h 4 ~ cS ) 

It follows that 



M < ft 



2-5 



< ft 3-c5 
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These errors only affect the fi equation at order h 4 ~ cS so fi is only affected at order 
h 3 ~ cS . Given this, the 9 equation is only affected at order h 3 ~ cS . Thus, the impact on 
9 is of size h 2 ~ c5 . In conclusion 



< h 2 - cS |/2| < h 3 - cS 



Thus 



it. 



Uz\\m < h 2 - c5 



Since u z in Theorem 1.1 in fact means u^, this completes the proof of Theorem 1.1 



Appendix A. Computations 



We shall carry out the computations of the ODEs appearing in (3.5) in Lemma 3.1 



and show that they are equivalent to (1.19), with errors of size 0(h ). This is carried 



out without making the symmetry assumption on the solution. When the even/odd 



symmetry assumption is imposed, we will carry out the integrals appearing in (1.19) 
and show that the ODEs claimed in the statement of Theorem 11.11 hold. 

Denote u z = 771+772- Let L = {1,2,5,6} denote the indices that refer to the left 
soliton and R = {3, 4, 7, 8} denote the indices that refer to the right soliton. The 
coefficient matrix of the symplectic form is 



[O'lr, 



A 







+ 0(h 2 



where the 0(h 2 ~) contributions come from a^ m with £ e L and m G R (and vice- versa, 
but of course ae m = —a m e). Fortunately, we do not need to compute these terms. 
Note that 

" I] 



A' 1 



+ 0(h 



2-\ 



—I 

In fact, we can substantially reduce the complexity of computation in applying 



Lemma 3.1 by observing that JH'{u z ) decomposes into terms parallel to M plus 
other terms which are 0(h 2 ~). To this end, we expand: 



H'(u z ) = H'( Vl ) + H'( m ) + #"(771)772 + #£(772)771 + 0(h 4 



where 



H p (u) 



\u\ 



Moreover, we have 



JH'(rf) = d v H(r])d a r] - d^H(T])d e T] . 



Hence, 



(A.l) 
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where 



b 2 = d^Hfa) , b A = d V2 H(r] 2 ) , b 5 = -d^H^) , b 7 = -d^H{r] 2 ) 



and all other bj = 0. Observe that (H"(r)i)r) 2 ,d z m z ) = 0{h A ~) for any t 6 R and 
(H p (r] 2 )rji,d z eu z ) = 0(h 4 ~) for any £ e L. Note further that for I e L (and hence 
d z m z = d z ir]\) we have 



(A.2) 



\H"(vi)V2,d z m z } = (r) 2 ,H"(r} 1 )d z tr} 1 ) = d z e{r} 2 , H'Jrjx)} 



Similarly, for i e R and (and hence d z iu z = d z eT] 2 ) we have 



(A.3) 



(Hp(r} 2 )r)i,d z j>u z ) = (771, #£(772) ^771) = d z t{r]i, H' p (rj 2 )} 



From (A.1),(A.2), and (A.3), we obtain 



d z eH{u z ) = -J2^^i + dAm,K(^)) + dAr j2 ,H' p (r ll )) 
i=i 



It follows that the equations (3.6) reduce to 



It suffices in this sum to discard 0(h 2 ) terms in a . Thus we obtain the equations 



' z 1 



i 2 



b' + d z ,{H'M 

b 2 + d z e(H , p (r Jl ) 

b 3 + d z r(H' p ( V2 ) 
b 4 + d z s(H' p ( V2 ) 

b 5 - d zl (H' p ( Vl ) 
z 6 = b 6 -d z ,(H' p ( Vl ) 
z 7 = b 7 -d z3 {H p {r] 2 ) 
z 8 = b 8 - dAH'M 



i 5 



ri2) + 0{h A - 

m) + o(h 4 ~ 

Vi) + 0(h 4 - 
Vi) + 0(h 4 - 

m) + o(h A ~ 
m) + o(h A - 
Vi) + o(h A - 
771) + 0(h A ~ 
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AH = 




+ d 6l (H'pirn ) 




+ 0(/i 4 ) 


di = 


+d vl H(r] 1 ) 


+ d vi (H p {rj x ) 


>m) 


+ 0[hr ) 






+ d 62 (H' p ( V2 ) 


,Vi) 


+ 0(h 4 -) 


d 2 = 


+d V2 H{V2) 


+ d V2 (H' p ^2) 


,Vi) 


+ 0(h 4 ~) 


9i = 






,m) 


+ 0(h 4 ~ 


Vl = 




- d ai (H'pfa) 


>m) 


+ 0(h 4 ~) 


9 2 = 




- d^(H p (r) 2 ) 


,Vh 


+ o(h 4 - 


v 2 = 




- 9 a2 (Hpfo) 




+ 0(h 4 -] 



We note that these equations hold in general, without assuming that the solution 
is even or odd. 

The next step is then to compute (H'(t] X ), r) 2 ) and (H' p (r) 2 ), rji). Let <f>(x) = sechx. 
We have 



(A.4) 



At this point we will make the even/odd assumption. In the even case, we may set 

(A. 5) 0, a, 9, v) = Oi, -ai, 9 X , -v x ) = (jjl 2 , a 2 , 6 2 , v 2 ) 

Then 6 X — 9 2 = . In the odd case, we may set 

(A. 6) 0, a, 9, v) = (fix, -a x , 9 X - vr, -v x ) = (p 2 , «2, 9 2 , v 2 ) 

Then 9 X - 9 2 = tt. 

In either the even or odd case, we find that (i — fi x — (i 2 — 0(h 3 ~), from which it 
follows that 

(A.7) // = A*i = A*2 = l + 0(h 2 -) 

Take a = in the even case and a = 1 in the odd case. We compute the equations for 
fi 2 , a 2 , 9 2 , v 2 by carrying out the appropriate derivative of (A.4), and then evaluating 
the resulting expression using (A.7), (A.5), (A.6). By residue calculus computations 
and asymptotic expansion, 



a(£,a 



dcf 



e -*X d) 3 (x-a) 



e -A»[4 + (2-4a)if + 



(x + a) dx 

_7r2 

~~6~ 



+ 2a-2a 2 )e + 0(e)}+0(e- 4a ) 
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/+oo 
e~ ix Z [0 3 ]' (x- a) <P(x + a)dx 
-oo 

= e- 2a [4+(6-4a)^]+0(/i 4 ~) 

We find that 

' f l = (-l) a Re[-ia] +0(h A -) 
a = + (-1) CT Re[+iaa + d ( a] + 0(h 4 ~) 

e = -fi 2 + -v 2 fi~ 2 + (-l) a Re[(iva + 3)a + v(d^a)} + Re(af3 - idtf) + 0(/i 4 ~) 

Z Z 

v = R e [-iva - (3] + 0{h 4 ') 

where a, d^a, and (3 are evaluated at £ = — 2v. 
Substituting, we obtain 

' ft = (-l) CT (8a - A)ve~ 2a + O(h^) 
a = + (-l) CT (-4a + \ix 2 )ve- 2a + 0(h 4 ~) 

= + i^" 2 + 18(-l) CT e- 2a + 0(h 4 -) 

Z Z 

v = _4(-l) CT e- 2a + 0{h 4 -) 

The system (/ta,i>) can be solved with error 0(h 2 ~); from which (a, i>) can be 
recovered with error 0(h 2 ~). At this accuracy the dynamics are comparable to 

{a — v 
v = -4(-l) a e- 2a 

Then \x can be solved with "explicit" order h 2 term coming from the order h 3 term 
in the equation for fi, and then 6 can be obtained with error of size h 2 . 
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